Network of edge states: random Josephson junction array description 
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We construct a generalization of the Chalker-Coddington network model to the case of fractional 
quantum Hall effect, which describes the tunneling between multiple chiral edges. We derive exact 
local and global duality symmetries of this model, and show that its infrared properties are identical 
to those of disordered planar Josephson junction array ( J JA) in a weak magnetic field, which implies 
the same universality class. The zero frequency Hall resistance of the system, which was expressed 
through exact correlators of the tunneling fields, is shown to be quantized both in the quantum Hall 
limit and in the limit of perfect Hall insulator. 



The QH transitions (!]] provide one of the best test- 
ing grounds for our understanding of quantum critical 
phenomena. The main advantages of this system are 
the availability of well-controlled samples, and the va- 
riety of field and gate-voltage-driven transitions within 
one sample. Unfortunately, it appears that current the- 
oretical understanding of the bulk QH effect is lagging 
behind the variety of experimentally observed features. 
Not only the transitions are found |^|,[| to have univer- 
sal scaling exponents and universal critical conductances 
HQ as suggested || basing on the bosonic Chern-Simons 
(CS) model ||, but their I-V traces at different filling 
fractions appear to be algebraically related in a surpris- 
ingly wide range of parameters. Most notable are the 
reflection symmetry J7J of I—V traces, and the quantized 
Hall insulator (QHI) phenomenon. 
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FIG. 1. Left: An illustration of the chiral edge network. 
Only the tunneling phases in saddle points (solid dots) remain 
in the effective tunneling model. Four-point measurement can 
be performed via the external gates G n coupled to long leads. 
Dashed lines outline some of the dual edges C„. Right: The 
simplest example of the finite-system duality between a QH 
droplet with a backscattering point and two separate droplets 
connected by a tunneling point. 

Theoretically, despite the attempted || rigorous analy- 
sis of the bulk CS model, the suggested universality || of 
the QH transitions has not been proven beyond the RPA. 
It was argued jlO| that both the charge-vortex duality 
and the QHI phenomenon (the semicircle law pl[ ) must 
be present simultaneously to account for the observed 
symmetry of I-V traces. Although the derivation of 
the semicircle law implies a duality between macroscopic 
current density and electric field, so far there was no 



quantum model to demonstrate both properties. Note, 
however, that the network of chiral edges, illustrated in 
Fig. ^a, shows both the correct scaling behavior (inte- 
ger regime, Chalker & Coddington and the QHI 
phenomenon (high temperatures, Auerbach & Shimshoni 
|l3|). Moreover, the effective action for a single tunnel- 
ing junction has an exact dual representation [^4|Jl5| . An 
important question is whether these features hold for the 
network of fractional edges at small temperatures. 

The goal of this letter is to construct the tunneling 
action for the network of chiral edge states. We derive 
local and global duality symmetries of the partition func- 
tion, and show that the infrared properties of the system 
are identical to those of the disordered JJA in a weak 
magnetic field, which implies the same universality class. 
We also express the four-point resistance of this system 
through exact correlators of the tunneling fields. 

Gapless edge excitations u = u(x, r) in the QH regime 
are described [fl6| by the imaginary-time action 



So = — I u 1 (iu + vu')dx dr, 



(1) 



in our units the edge wave velocity v — 1. The action 
remains quadratic |ll],|l7j in the presence of Coulomb 
forces; the interaction is introduced by the tunneling 



St 



i<j<2N 



dr 53 A l cos(gk(j) t3 ), 
k=l 



(2) 



where g — *Jv for the quasiparticles' tunneling between 
the points Xi and Xj through the condensate with the 



filling fraction v — (2m +1) , or g— >g = X/ \fv for tun- 
neling of electrons through an insulating region. The tun- 
neling amplitudes , set by the details (T^j of the self- 
consistent potential near saddle points a=ij, l<a<N 
are treated as phenomenological parameters. Although 
the argument of the cosine in Eq. ^| is usually written 
as the difference <j>ij = Uj — Ui, we must remember that 
tunneling connects two separate points, and there may 
be an additional gauge field contribution |L8| . 

The non-linear action (||) depends on values itj of field 
u in a discrete set of points xf, its fluctuations in all other 
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points can be integrated out. Leaving the arguments <p a 
of the tunneling terms as the only independent variables 
(cf. gauge- invariant phases for Josephson junctions), we 
write the most general form of the tunneling action 



S[(/> a ] — T ^2 j-<t>c,Kap(t>i3 + S t , 

u=2nnT 



(3) 



with Matsubara frequencies u> and 4>a{^) = 4>a.(—u). Dif- 
ferent tunneling points are linked by the elements of the 
coupling matrix K a ^(uj) = —Kp a (—oS), determined by 
the geometry of edge channels; the factor lu/Att was sep- 
arated for notational convenience. 

Before specifying the form of the coupling matrix K, 
let us discuss general properties of the action (^): 
(i) In the absence of any tunneling the matrix K is 
uniquely determined by correlators of the fields cf> a , 



(4) 



Therefore, the matrix elements between the tunneling 
points located at otherwise disconnected edges vanish, 
(ii) In the limit of infinite frequencies, the coherence be- 
tween different tunneling points is also lost, while the 
matrix K acquires its asymptotic form 



K a p(uj) 5 a psgn{ui) 



(5) 



(iii) The action (|J) describes a complicated system of 
non-linear damped anharmonic oscillators. Since their 
mutual coupling is linear, each of these oscillators may 
be regarded as a particle 



& = |j:I>[tfuMI& 



in its own periodic potential and external field Vi(u>) = 
^2 a ^ 1 Ki a {ijj)(j) a {uj) due to all other tunneling points, 
whose dynamics is determined by the remaining terms 
S'[(f>2, • ■ ■ , 4>n\ of the original action. The partition func- 
tion can be calculated in two steps, 

Z= Je- s V[cj)}= Jvfo...V<j) N e~ s ' [ ^-^ N] Zi [Vi] , 

where Z\ [Vi] = JVfa exp (—Si). This one-point parti- 
tion function can be rewritten p3 in terms of the dual 
variable 9\ with the action 



K n (u>) 



|0i(w)+Vi(w)| 2 + $[0i(t)], (6) 



where the tunneling part St[9] has exactly the form (|^) 
with modified tunneling coefficients and additional 
replacement g — > g = g~ 1 . This duality represents |]l5f a 
freedom to describe the same junction in terms of weak 
tunneling or strong backscattering as illustrated in Fig. [I]. 
If applicable to larger systems, such duality must severely 
restrict the form of the matrix K. 



We illustrate these properties by calculating the cou- 
pling matrix for systems in Fig. |l|b,c. For the single-edge 
geometry in Fig. |l|b, the total charge Q = Jp(x) dx — 
g (ul—uq)/2it is conserved, and it is sufficient to consider 
periodic configurations u x +l = u x , where L = L\ + L 2 is 
the total edge length. The integration over fluctuations 
of the Gaussian field u is performed by evaluating the 
action ([j]) along the classical trajectory connecting the 



points Uj 



S [uj,Uj-i 



and Uj-i, j 



1,2, with the result 

(ujSj - uj-i) 
J— I ■ 



(7) 



Here the distance factor Sj = exp(ujLj) depends on the 



distance along the edge Lj = Xj 
of the coupling matrix K is 



Cj-i 



The only element 



K 



(ib) 



s\s 2 - 1 



j(L 1 +L 2 ) _ ^ 



(s x - l)(s 2 - 1) ~ (e" L i-I)(e^_i)- ( 8 ) 

Clearly, the property (ii) is satisfied: at high enough 
frequencies, |w|£i,2 ^> 1, the expression (||) reduces to 
sign(w), up to exponentially small corrections. 

The two-edge geometry shown in Fig. [j]c is homomor- 
phic to that of QH liquid on the surface of a cylinder 
investigated by Wen |l6j , who emphasized that mere col- 
lection of independent edges does not give its complete 
description |2(J . This is related to the fact that tunneling 
destroys the charge conservation at individual edges, and 
the periodic boundary conditions arc no longer valid. It 
turns out that the free boundary conditions (ul indepen- 
dent of Uq) combined with the "boundary terms" 

S C nd = T/(8tt)J2 lu ^{ul +uq){u l -U ), (9) 

which must be added to the action (|l]) for every separate 
edge, restore the correct Hilbert space. Particularly, for 
the system in Fig. 0c, after integrating out all Gaussian 
modes, this prescription leads to the effective action 

uj\9u,\ 2 (e wil -l)(e wL2 -l) 



S = TJ2- 



-in 



s t 



(10) 



which satisfies the finite-size duality requirement (iii). 

Although the same prescription works for larger net- 
works, in systems with more than four tunneling points 
direct calculation of K becomes too bulky. Instead we 
use Eq. (^) , which implies that the element of the inverse 
coupling matrix K~p is independent of the presence of 
the tunneling points 7 ^ a, (3. The diagonal elements of 
matrix K~ x can be always deduced from Eqns. (||) or jlfj| ) 
by substituting the appropriate lengths for Li and L 2 , 
while the off-diagonal ones can be derived by construct- 
ing effective actions with only two tunneling points. 




1 ~ 3 TV 

FIG. 2. Non-trivial geometries with two tunneling points. 
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The coupling matrices for two-tunneling-point geome- 
tries shown in Fig. are given by the expressions 



£(2a) _ 



S1S2S4 — 1 
(si -l)(s 2 s 4 -l) 

S4 

— 1 



32^3^4 — 1 

(*3— l)(aa»4— 1) 



(S 1 S2-1)(S 3 S 4 -1) 

23 1 2L. 



1 — S3S4 1— S1S2 

3lS 2 S3S4 — 1 



1 — S3S4 1 — sis 2 (sis 2 — l)(s 3 s 4 — 1) 



(11) 

(12) 



where s< = exp(wLi). All other geometries with two tun- 
neling points can be generated from these by replacing 
one or both tunneling junctions by their dual(s) as il- 
lustrated in Fig. [l]b, c. Our calculation shows that the 
coupling matrices for thus obtained configurations are 
immediately related to Eqns (|TT|), ( |T2] ) as prescribed by 
Eq. (^): the modified coupling matrix is given by the 
coefficients of the Legendre-transformed bilinear form 



4tt 



(13) 



This implies that under the global duality transforma- 
tion, which replaces all tunneling junctions by their du- 
als, the coupling matrix is inverted, K = K . Combin- 
ing this statement with Eq. (||) , we find that the original 
matrix K is given by the correlators of the dual fields 



/ ~ \ 

K a(3 = \SaOp)~ x = 



A=0 



(14) 



in the limit of perfect transmission (zero backscattering) 
in every tunneling point. 

Although the dual form of the model is the result of 
a non-local change of variables in the path integral, the 
tiieafc-coupling limits of the original model and its dual 
version can be interpreted as describing the same system 
at mutually dual filling fractions, with the total areas 
of the condensate and the depleted regions interchanged. 
This interpretation is correct as long as the probability 
distribution of the bare parameters of the original model, 
and the model at the dual filling fraction expressed in 
dual variables are similar, which is expected if the distri- 
bution of the disorder potential is symmetric. Under such 
duality "in average" g—^g^ 1 , and therefore at g ^ 1 this 
transformation is not an exact symmetry of the prob- 
lem. This is obvious in the limit of large frequencies, 
where the off-diagonal elements of the coupling matrix 
K vanish, and the scaling dimensions of the tunneling 



amplitudes coincide with their values £k = 1 — k 2 g 
for isolated tunneling junctions. 

At small frequencies, however, the interference be- 
comes important, and the scaling behavior crosses over 
to mesoscopic regime where the quantization of the drift 
orbits is relevant. It turns out that the quadratic part of 
the action in this limit can be written as a simple sum 

u: — >U — Lip 

C C 



2„2 



where = ^acc^" IS ^ Qe directed sum of the tunnel- 
ing phases <j) a along the dual edge C of length created 
by the global duality transformation. The presence of 
such terms in the action imply that variables can have 
only finite r.m.s. deviation, which ensures that the aver- 
ages = ^2 aC c^ a ) ~ vanish. These constrains 
can be resolved by writing the phases <j) a as the gradi- 
ents (j) a = gA a V + A a of the local potential V = Vc, 
associated with every independent edge C, up to some 
time- independent phases A a . In this representation the 
model becomes identical to a disordered JJA in external 
magnetic field. Similarly, in the dual representation of 
the model the sum 0c = J2acc ® a can ^ e associated with 
the total current of composite bosons entering the en- 
closed area; this current also vanishes at zero frequency. 

The performed analysis of the infrared limit reveals 
a deep analogy between the CB in disordered QH sys- 
tems and bosons near disorder-driven superconductor- 
insulator (SC-I) transition. After the system is separated 
into weakly coupled phase-coherent areas, the terms re- 
maining in small-frequency expansion of the action de- 
pend on lengths L a in a highly non-universal fashion. 
This is used to absorb the Luttinger liquid coupling con- 
stant g, which renders the tunneling model independent 
of the filling fraction u, leading to the conclusion that all 
QH transitions must be in the same universality class, 
associated with disordered SC-I transition. 

This universality follows only from the mapping be- 
tween the partition functions, and it does not imply the 
identical transport properties. The chiral nature of the 
current-carrying states in the edge network model is re- 
sponsible for its large Hall resistance. This can be readily 
seen in the integer regime using the Biittikcr-Landauer 
pl[ formula relating the four-point resistance tensor with 
the matrix T nn i — \t nn >\ , n,n' — 1,...,4 of transmis- 
sion probabilities between the incoming and outgoing 
channels as illustrated in Fig. |l|a. Although for an ar- 
bitrary unitary scattering matrix t nn i the Hall resistance 
i?H = {Rxy — Ryx)/2 is not necessarily quantized, in the 
classical regime, which is characterized by the absence 
of the quantum interference between different paths con- 
tributing to t nn i , i?H = ±1) in agreement with Ref. p3| . 

Within the proposed model, the infinite ideal leads 
C n can be coupled with the potential gates G n using 
the action Sc n — 2tt i g^ 1 jdr <d n W n , where the phase 
0„ = X)qcC @<x ^ s proportional to the total charge trans- 
ferred to the lead C n , and W n is the potential of the 
corresponding gate G n . In the linear responce regime 
the current to the lead n is given by the expression 



Z„(a;) = ujg- 2 {Q n Q nl )W n ,{Lo), 



(15) 



while the equilibrium potential differences between the 
opposite leads can be evaluated by decomposing them 
into individual chiral channels, e.g. 

2g{V 1 -V 3 )+g-\l2-I i )=Ii-l2-i 3 +ii- (16) 
The currents 
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are denned at the external edges C n (dashed lines in 
Fig. in the dual representation of the model. The 



r.h.s. of Eq. (17) can be related to the potential differ- 
ence between the phase-coherent regions of leads 1 and 3 
in the effective JJA model, which leads to the formula 



R=-g~ 2 ^+R h , 



for the resistance tensor as the sum of the quantized Hall 
resistance and the part Rb associated with the trans- 
port of bosons in the JJA. An analogous decomposition 
P — Pquant + Pb wa s previously derived H for resistivities 
in the RPA approximation; we believe this expression is 
exact here because the charges are essentially localized by 
disorder and the fluctuations of the Chern-Simons field 
are suppressed. The usual relationship pb — > implied 
by the particle-vortex duality |2^,|^,^] in bulk models, is 
valid in the considered model only after averaging over 
the disorder. 

Eliminating auxiliary potentials W n from Eqns. (|l5|) 
( |l7| ) and an expression similar to Eq. ( 16) for leads 2 and 
4, we finally obtain 



(17) classical resistor network, which demonstrates jl3| the 
QHI phenomenon, in general this behavior may not be 
present. Finally, since the constructed model involves 
no additional approximation compared with the full chi- 
ral Luttinger network model, it can be used to simulate 
tunneling experiments in the fractional QH regime. 

It is our pleasure to acknowledge beneficial discussions 
with S. Chakravarty, E. Fradkin, I. Gruzberg, C. Kane, 
S. Kivelson, D.-H. Lee, C. Marcus, E. Shimshoni and S.- 
C. Zhang. 



(18) 



i?b : 



: . 9 2 (i+s)[iL»($e) (ee)- 



where the matrices 



S = 



1 






-1 



D 



-1 
1 



(19) 



(20) 



represent the configuration of currents and potential dif- 
ferences in the four-lead measurement. Note that after 
the rescaling g $ — > $, Q/g — > 6 used to construct the 
effective JJA model, the expression ([l9]) loses its explicit 
dependence on the filling fraction v = g 2 . 

The correlator ($6) in Eq. (jl?]) is the linear combina- 
tion of the expressions 



OaV/3) 



2?A 



a/3" 



7/3 = 



7-yVp) 



According to Eqns. (Q) and ([Tj), these vanish identically 
if the tunneling is absent in either the original (A Q = 0) 
or the dual (A Q = 0) representations, which leaves only 
the quantized part of the Hall resistance. 

To conclude, we constructed the effective tunneling 
model which generalizes the Chalker-Coddington model 
|L2| to the fractional QH effect. The partition function 
of this model can be approximately mapped to that of 
the disordered JJA, which implies the universality of the 
quantum Hall transitions. The model always has an ex- 
act dual representation, but it maps to the same system 
at different filling fraction only in the limit of small fre- 
quences and after averaging over disorder. The associ- 
ated relationships between the filling fractions and the 
transport coefficients are identical to those obtained in 
the bulk CS model. Although in the limit of large tem- 
peratures the constructed model becomes equivalent to 



[1] R. Prange and S. M. Girvin (ed.), The Quantum Hall 

Effect, (Springer- Verlag, NY, 1990). 
[2] H. P. Wei, D. C. Tsui, M. A. Paalanen, and A. M. M. 

Pmisken, Phys. Rev. Lett. 61, 1294 (1988); L. Engel, 

H. P. Wei, D. C. Tsui, and M. Shayegan, Surf. Sci. 229, 
13 (1990); H. P. Wei, S. W. Hwang, D. C. Tsui, and 
A. M. M. Pruisken, ibid. 229, 34 (1990); S. Koch, R. J. 
Haug, K. von Klitzing, and K. Ploog, Phys. Rev. Lett. 
67, 883 (1991); L. Engel, D. Shahar, C. Kurdak, and D. 
Tsui, ibid. 71, 2638 (1993). 

[3] L. Wong, H. Jiang, N. Trivedi, and E. Palm, Phys. Rev. 

B 51, 18033 (1995). 
[4] D. Shahar et al, Phys. Rev. Lett. 74, 4511 (1995). 
[5] S. Kivelson, D.-H. Lee, and S.-C. Zhang, Phys. Rev. B 

46, 2223 (1992). 
[6] S. M. Girvin and A. H. MacDonald, Phys. Rev. Lett. 

58, 1252 (1987); S. C. Zhang, T. H. Hansson, and S. 

Kivelson, ibid. 62, 82 (1989). 
[7] D. Shahar et al, Science 274, 589 (1996). 
[8] D. Shahar et al, Sol. Stat. Commun. 102, 817 (1997). 
[9] L. Pryadko and S.-C. Zhang, Phys. Rev. B 54, 4953 

(1996) ; E. Fradkin and S. Kivelson, Nucl. Phys. B 474, 
543 (1996); L. Pryadko, Phys. Rev. B 56, 6810 (1997). 

[10] E. Shimshoni, S. L. Sondhi, and D. Shahar, Phys. Rev. 

B 55, 13730 (1997). 
[11] A. Dykhne and I. Ruzin, Phys. Rev. B 50, 2369 (1994); 

I. Ruzin and S. Feng, Phys. Rev. Lett. 74, 154 (1995). 
[12] J. Chalker and P. Coddington, J. Phys. 21, 2665 (1988). 
[13] E. Shimshoni and A. Auerbach, Phys. Rev. B 55, 9817 

(1997) . 

[14] A. Schmid, Phys. Rev. Lett. 51, 1506 (1983); F. Guinea, 
V. Hakim, and A. Muramatsu, ibid. 54, 263 (1985). 

[15] P. Fendley, A. Ludwig, and H. Saleur, Phys. Rev. B 52, 
8934 (1995). 

[16] X.-G. Wen, Phys. Rev. B 41, 12838 (1990); X.-G. Wen, 
Int. J. Mod. Phys. B 6, 1711 (1992). 

[17] H. Yi and C. Kane, Phys. Rev. B 53, 12956 (1996). 

[18] This contribution was not taken into account in Ref. jl^] . 

[19] J. Jain and S. Kivelson, Phys. Rev. B 37, 4111 (1988). 

[20] Otherwise one gets unphysical selection rules for tunnel- 
ing through a fractionally charged antidot [V. Pokrovsky 
and L. Pryadko, Phys. Rev. Lett. 72, 124 (1994).] 

[21] M. Biittiker, Phys. Rev. Lett. 57, 1761 (1986). 

[22] M. P. A. Fisher, Phys. Rev. Lett. 65, 923 (1990) 



